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I. INTRODUCTION 

The mechanism of unconventional superconductiv- 
ity remains a major source of interest for theo- 
rists. Anisotropic crystals, in particular organic 
ones, are expected, both from a theoretical point 
of view and experimentally, to show unconventional 
behaviour. Among the many materials, which 
are currently studied, SrjA^CaxCu2AOii^ Q or 
{LaY)y{SrCe)ii-yCu2iOii^ [j| show very anisotropic 
structures, and can be well represented by a single lad- 
der. Indeed, unconventional superconductivity has been 
observed in these materials 0, [1]; in particular, super- 
conducting (SC) phases are found in the vicinity of an- 
tiferromagnetic ones, sometimes with a gapped phase in 
between^?]. These observations indicate the possibility of 
coexistence of magnetism and superconductivity. 

The materials we have studied in this work belong 
to the class of strongly correlated quasi-one-dimensional 
systems of electrons. One important peculiarity of these 
systems is the existence of nearly degenerate phases with 
different symmetries and, therefore, the possibility of 
very rich phase diagrams, commonly observed experi- 
mentally and often predicted theoretically. In quasi-one- 
dimensional systems, it is well established that two kinds 
of Fermi surface instabilities can occur: (i) Cooper insta- 
bilities, including singlet or triplet superconducting ones; 
(ii) Peierls instabilities, including Spin Density Wave 
(SDW) or Charge Density Wave (CDW) ones. A very 
important and quite general property is that these two 
kinds of instability are coupled, because of the topology 
of the phase space in one dimension. The questions about 
how these instabilities compete and sometimes coexist 
are very important for the physics of these systems. Be- 
cause of this instability coupling, such a discussion can- 
not be done in a Fermi liquid approach, but requires more 
sophisticated methods. 

Different approaches have been used, to understand 
unconventional superconductivity, like Monte Carlo sim- 
ulations with the t—J model[^, Q or with the Hub- 
bard modelfiol. [Tl|. or exact diagonalizatio nf]3 . DMRG 
method[l^ or variational approach p^. etc.jl5l|. 

An important step forward has been taken with the use 
of the Renormalization Group (RG)[T6|. Not only can 



these calculations predict a SC phase [TtI, [l^ . but they 
give a new interpretation of this unconventional mecha- 
nism: it results from the competition between the Cooper 
channel (formation of pairs of electrons) and the Peierls 
channel (formation of of electron- hole pairs) [l^ . 

The RG method is a fixed point method, its applica- 
tion in condensed matter has a severe drawback: the RG 
flow is always diverging, so that no exact fixed point can 
be obtained; in other words, it is impossible to calculate 
the renormalized parameters of these systems. Never- 
theless, one can calculate the phase diagram, by exam- 
ining which susceptibilities are diverging (i.e. are unsta- 
ble) and which are remaining finite: the processes corre- 
sponding to non divergent susceptibilities are negligible 
compared to those corresponding to divergent ones. 

This paper is devoted to the study of a ladder sys- 
tem, which consists of two coupled chains of atoms, the 
intrachain coupling is written the interchain one 
with <C t\\. We use the Hubbard model, which has 
been widely studied by theorists [20|. though its c omp lete 
analysis hasn't yet been achieved. M. Fabrizio[21] has 
previously calculated the phase diagram of the ladder by 
a two-loop expansion using RG equations. He obtains 
a very rich diagram, with Hubbard parameter U rang- 
ing over [0, IS-KVf] (vf is the Fermi velocity) and t± over 
[0, 1.2Ao] (Ao is the half band width), though its vahdity 
is somehow questionable, since U /vf is the parameter of 
this expansion. If one focuses on range U G [0, 27rvf] , Fab- 
rizio predicts a superconducting phase, which he named 
phase I; in this phase, the RG flow of susceptibilities 
shows several divergences: the SDW channel coexists 
with the superconducting one. 

We proved recently, for small values of the interchain 
interaction t±, the existence of an extra SDW phase, in 
this region of parameters, by including ify dependence 
of the couplings |22]. This phase is characterized by the 
flow of all superconducting susceptibilities, which remain 
flnite, while SDW ones diverge. These calculations have 
been performed with a fixed Fermi surface. This work 
also established the importance of high energy processes 
(like the backward intcrband scattering g^, see below) 
during the RG flow: although these processes die before 
the flow becomes divergent, they prove eventually influ- 
ential. 
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In this paper, we will discuss the effect of the renormal- 
ization of the Fermi surface, in the line of these i^n depen- 
dant RG calculations. One of the questions is whether 
our results, in particular the existence of a SDW phase, 
are valid or not. The answer is fortunately yes. 

In the last decade, RG methods have achieved very so- 
phisticated schemes: here, we use either the One Particle 
Irreducible (OPI) scheme, following; Ref. [H, H^li or the 
Wick-ordered one, following Ref. [23, and calculated the 
scatterings in a one-loop expansion. The renormalization 
of the Fermi surface remains valid, in this approximation. 

A remarkable result is that the phase diagram becomes 
scheme-dependant. This question was first addressed by 
H. Schulz, who argued that high energy processes would 
be influent in specific cases: this implies that the way 
they are included in the RG calculation would matter 
The response it receives here contradicts the usual opin- 
ion, shared by a number of specialists, that all schemes 
are equivalent and give identical results. 

We will first describe the model (section [ll]) and the 
RG equations (section [1111) . then discuss the choice of the 
RG scheme (section [TV|) and analyse our results (section 



II. MODEL 

In a ladder, there are two separated bands in the dis- 
persion diagram (0: binding and tt: antibinding) , be- 
cause of the Coulombian interaction between the chains. 
In other words, in the K± direction, there are only two 
physical points, O and ir/b (b is the interchain distance). 
There are four Fermi points {—kfo, —kf^^, kfT^, kfo) in 
the i^ii direction (see Fig. [T|). We will simply note K, for 
the momenta in the direction (and k will always be 
the relative momentum to a given Fermi point). 

The Fermi surface gap is defined as Akf = kfo ^ kf-j^. 
From Luttinger theorem, kfo + kf^^ is constant, so Akf 
is the only Fermi surface parameter. It relates t± the 
interchain interaction by Akf = 2t±/vf. 

The kinetic Hamiltonian is linearized around the Fermi 
points 26] with a single Fermi velocity Vf, and writes {R 
reads right moving particle and L left moving one): 

Hoin = E ( E(^^ - kfo)Rl{K)R^AK) + 

a- K 

{K - kfMAK)R.AK) + {K + kfo)Ll^{K)Loa{K) 
+ {K + kf^)Ll,{K)L,,{K)] 
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Figure 2: Schematic definitions of the couplings 9 



The interaction Hamiltonian writes 
1 

N 



Ki,K2,K'j,K^ 0-l,CT2 
Kl+K2=K'j+K^ 



-^K2Cr2^Kicri + S4-t'K'^cri^K^cr2^K2Cr2^Kicri 

-fgi(Ki,K2,K;„Kl)<,^^Lt,,^^i?j,^^^Lj,^^^ 
+S2(Ki,K2,K^,Kl)i?t^,^^Lt,,^^Lj,^^^i?j,^^^ 

in which Sq, is the two-particle coupling, and we have 
used the g-ology representation. More precisely, there 
are 8 different couphngs go, ffTi-, ff/o, 5/7r, ffto, QtT,: gbo and 
gbTT, corresponding to the interaction processes shown in 
Fig. [2] (in this way, all K± dependence of the couplings 
is included in the symbolic names, whereas all i^n de- 
pendence is given in their arguments, see more details 
in Ref. [221]), plus the Si = SiRRRR) = 9{LLLL) cou- 
plings which are not renormalized in a one-loop expan- 
sion. At the beginning of the RG flow (A = Aq), all 
scatterings Sq. are set to U, the Hubbard constant, thus 
one simply gets S4 = U. 



III. RG EQUATIONS 

The RG equations for scatterings go, g-ir, g/o, gb-n 
express their derivative as the sum of two terms: the first 
term is usually called Cooper term, since it comes from an 
electron-electron diagram; the second one is called Peierls 
term, since it comes from an electron-hole diagram. One 
can write, in a generic way. 
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Figure 1: The 2-band dispersion in || direction 



E 

K',K" 



g(Ki, K^, K^j, k;)s(k'/, K^', K^', + 



3 



K' 



K ; K 



Figure 3: One-loop tadpole diagram 
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S(K;, Ki, K^, K^)g(K'/, K.I K.I K.'i) (1) 



in which C and T are coefficients (C stands for Cooper 
term while T stands for Peierls term); explicit and de- 
tailed sums are given in Appendix B of Ref. (2^ ; £ is the 
flow parameter (the half band width is A = Aq e~^). 

For all couplings, except gho and gh-K, we get, in the OPI 
scheme, 6 = 1/ {A + 2\ki + K2\) and, in the Wick-ordered 
one, e = l/{A-2\ki+k2\) (here K = VfK/A); we get, in 
the OPI scheme, ? = l/(4+2|^i-X3|) and, in the Wick- 
ordered one, y = 1/(4 — 21^1 — K^l). In fact, the generic 
expression ([!]) does not apply to couphngs gt: the Cooper 
term splits into two terms, one with the same C factor, 
one with a special factor C^p = 1 / (4 + 2 1 + ^2 ± 2 Afc/ 1 ) 
for the OPI scheme and &p = l/(4-2|Xi + X2±2Afc/|) 
for the Wick-ordered one {Akf = VfAkf/A); the Peierls 
term splits into two terms, one with the same T factor, 
one with a special factor T^'p = l/(4-t-2|^i-7^3±2Afc/|) 
for the OPI scheme and TP = l/(4-2|^i -i^3±2Afc/|) 
for the Wick-ordered one (± reads + for gt,o and — for 

The RG equation for Ak f is obtained through the two- 
loop expansion of the self-energy S, following a standard 
calculationjUi^,!!!!. Let Go = Z/{-iuj + Vf{K - kfg + 
/x)) be the free Right propagator of the band 9 {9 — 0,n), 
and fi the chemical potential, one can write 
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G„ 



^R0 = 5G-^ = -{Svf{K-kfg)-VfSkfe + dn) ^dZ . 

In Fig.[3l we show the tadpole diagram, corresponding 
to a one-loop contribution in this expansion of S: af- 
ter all simplifications (one has to subtract carefully the 
contribution of Sfj,), one gets, in this one-loop expansion. 



57r2 - .902 + 5/71-2 — 5/02 



<)iTl — 501 + 5/7rl - 5/01 



It is obvious, in this formula, that Akf depends on K^^, 
however, this dependency gives very small variations and 
can be neglected. 

There are two different two-loop diagrams, represented 
in Fig.m The first one (a) gives no contribution, and the 
second one (b) (sunrise) gives three. 



Figure 4: Two-loop tadpole diagram 



Two of them give logarithmic terms; in fact, these two 
contributions can be deduced one from the other using 
symmetry A A' (see notations in Ref. [2^). 

The only contribution containing S4 is 
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We have calculated the complete expression of the two- 
loop expansion, including i^n dependence, both in the 
OPI scheme and in the Wick-ordered one. Except for 
the last Si contribution, these expressions depend on the 
RG scheme. It would be too fastidious to explicit all 
contributions: instead, let us skip all K^^ dependence. 
We get the following contribution: 



(5(Afc/)J 2-loop 
except 3 4 

-2Ao- 



iZSA / Vfk-A Vfk + A 
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Vfk 
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for each 50, 5/o, 5to (for which ± reads +), g-^, gfn, and 
gtTT (for which ± reads — ). This is similar to previous 
calculations (Ref. [2^ for the OPI scheme, Ref. [s^l for 
the Wick-ordered one), but we would like to emphasize 
one major novelty: for g^Q and gbn, the second factor is 
modified and writes 

Vf{k±2Akf)-A Vf{k±2Akf) + A 



\vf{k±2Akf) - A\ \vf{k±2Akf) + A\ 

-2Ao-vf{kT'2Akf) 2Ao~Vf{kT'2Akf) 
°^^-A-\A-Vf{kT2Akf)\ A+\A + Vf{k±2Akf)\ 

in which the ± reads as in the first factor. 
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To end with technical details, let us explain the approx- 
imations used in the RG equations. First, all scattering 
S depend on three arguments (^1,^2,^3), which are re- 
placed by their 2piAkf {pi G Z) approximation. This 
is generalized to all other couplings. Second, the list of 
all functional couplings S((fcij ^2, fcs) is truncated by set- 
ting \pi\ = 2, 3 or 4. Extra couplings are replaced by 
the closer element in the list using symmetry preserving 
relations (cf. Ref. Last, couplings S((fci, fe, ^3) in 

which some |fci| ^ 2/S.kf are replaced by U , the Hubbard 
constant (this happens when the initial half band width 
Aq S> A/c/, i.e. for small values of it mostly arises 
from the logarithmic contributions). 



IV. CHOICE OF THE RG SCHEME 

It is not the place here to derive the RG equations for 
the OPI scheme [3l|, nor for the Wick-ordered one [3^. 
What matters here is that one can express the RG flow 
in terms of couplings S, Fermi gap Afc/, Fermi velocity Vf 
and renormalization factor Z . As far as we will not dis- 
tinguish VfQ and i;/7rj we need not discuss the renormal- 
ization of Vf and Z, which only induces a global scaling 
of the other couplings, subsequently we will forget these 
parameters. 

To get the RG equations, one expands diagrammati- 
cally all couplings, as in the Cauchy expansion in U /vf. 
For a given energy scale A, one of the inner energies is 
integrated in the range [A — 5A, A -I- (5A], where JA is in- 
finitesimal; in the OPI scheme, all other inner energies 
are integrated over [A, Aq]; in the Wick-ordered scheme, 
they are integrated over [0, A] [33] . OPI and Wick-ordered 
schemes not only differ according to these rules, they also 
give different C and f factors, as explained before. 

From a theoretical point of view, both schemes should 
converge to the same fixed point, however, the RG flows 
are divergent and therefore never reach the fixed point: 
the integration of energy is incomplete; therefore, it is 
crucial to choose whether one will integrate over UV en- 
ergies first (i.e. \E\ > A, as in the OPI scheme) or over 
IR energies first (i.e. \E\ < A, as in the Wick-ordered 
one)!!!. 

This choice is expected to be more influential when 
high energy processes are taken into account. Within 
the Wick-ordered scheme, such processes participate in 
the RG flow at the very beginning, when A = Aq, but 
they are skipped when £ is increased. In the OPI scheme, 
they are always taken into account. 

In the ladder system, there is one such process, corre- 
sponding to the backward interband scattering gi,. It is 
indeed a high energy process, only permitted for \E\ > 
2/S.kf. Within the Wick-ordered scheme, this contribu- 
tion is suppressed for £ > ln(2w/Afc//Ao). After these 
considerations, one could expect that the RG calcula- 
tions performed with a fixed Fermi surface would bring 
different results, depending on which scheme is chosen. 
However, the weight of the gi, contribution is proportional 



to e, e'^P, T or y^P. In the OPI scheme, all these terms 
vanish as 1/(1 + —x^), when £ is increased, so gb mostly 
contributes to the RG flow at the beginning, as in the 
Wick-ordered scheme. We have indeed performed both 
calculations and found a difference which is meaningless 
and negligible (35}. 

However, if the Fermi surface is correctly renormalized 
during the RG flow, in the case when Akf — > as far 
as A — > 0, the weight of the gh contribution keeps finite 
during the RG flow, in the OPI scheme, whereas it is 
still suppressed for large values of £ in the Wick-ordered 
one; so, the results of RG calculations should prove sig- 
nificantly different, using one or the other scheme. 

In our opinion, the choice of the OPI scheme is more 
convenient, because, in the Wick-ordered scheme, the 
weight of high energy processes is underestimated. This 
is, in particular, the conclusion of C. Nickel, who has per- 
formed a careful comparison of different RG schemes (see 
subsection 3.4 of Ref. H). 

There is another indication that it is more correct to 
use the OPI scheme: in his pioneer work with H. Schulz, 
D. Zanchi[36j| has proved that some terms in the 3-loop 
expansion induce an integration of energies \E\ > A. 
Therefore, the OPI seems the only self-coherent scheme, 
when one tries to include further terms in the perturba- 
tive expansion. 



V. RESULTS 

The results of these calculations confirm those of 
Ref. [131, done with a fixed Fermi surface (i.e. Akf 
was kept constant). We find two distinct regions; in 
the SDW region, no superconducting susceptibility is di- 
verging, while SDW ones are (see Fig. \5\ (a)); in the SC 
region, both are diverging, but the superconducting sus- 
ceptibility always dominates (see Fig. [5] (b)). However, 
the intermediate region, described in Ref. in which 
superconducting susceptibilities, although diverging, are 
not dominating, vanishes completely. 

In Fig. [51 the two phase diagrams are presented, ac- 
cording to the choice of the RG scheme. Let us repeat 
that, when AA:^ is not renormalized, both schemes give 
almost the same phase diagram[35]. Here, on the con- 
trary, one observes that the SDW region quantitatively 
depends on the RG scheme. Indeed, for all values of U, 
except very small ones, in the Wick-ordered calculation, 
a constant critical value t±c can be defined, which sepa- 
rates the SDW and the SC regions. In the OPI calcula- 
tions, the evolution of this critical value tj_c is smoother, 
with a linear part of slope 8.6 at small U. On the 
whole, the critical line t±ciU) which separates both re- 
gions differs quantitatively, except for small values of 
C/HJ. With the Wick-ordered RG scheme, the SDW 
area is reduced by a factor 3, compared to the result of 
the RG with the OPI scheme. 

The difference of results coming from the choice of the 
RG scheme has already been suggested by several au- 
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(a) 



(b) 



Figure 5; Flow of the susceptibilities for 2t±/Ao — 1.4 (a) or 
2tx/Ao = 1.5 (b), and = 0.5. xf*^''''(0) is the intraband 
singular SC susceptibility of d symmetry, X^'^^^H^) is the in- 
traband triplet SC susceptibility of / symmetry, Xc^(O) is 
the intraband CDW susceptibility, Xs^(O) is the intraband 
SDW susceptibility, Xc^if) is the interband CDW suscepti- 
bility and Xs^(f ) is the interband SDW susceptibility (the 
difference between site and bond susceptibilities as well as the 
symmetry classification are explained in part III of Ref. [13]). 
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Figure 6: Phase diagram, versus parameters t± and U — 
the central area belongs to the SDW phase, accord- 
ing to the OPI scheme, and to the SC one, according to 
the Wick-ordered scheme. The parameter n indicates range 
[— 2nAfc/, 2nAfe/] in which K\\ dependency is exactly taken 
into account. We did not distinguish the curves for n — 2 
and n = 4 in the Wick-ordered scheme, because they hardly 
separate de visu. 



thors. It was, in particular mentioned in Ref. ^2^. 

The behaviour of the Fermi surface, during the RG 
flow, brings no surprise. Let us first present the results 
in the SC phase, then in the SDW one. 

In the SC phase, Afc/ increases slowly, while not di- 
verging. The flow diverges at some Ac, and Afc/c is the 
final value of Afc/. The numerical values of Akfc/Ao are 
not realistic, however the general trend is very satisfac- 
tory and indicates that the binding/ antibinding separa- 
tion is necessary to the existence of superconductivity. 
We believe that, if one would include a greater number 
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Figure 7: Flow of the Afc/: in the SDW phase, (a) using OPI 
RG calculations, or (b) using Wick-ordered RG calculations; 
in the SC phase, (c) using OPI RG calculations, or (d) using 
Wick-ordered RG calculations. 



of chains in the model, one would obtain more realistic 
values for Akfc- 

Let us emphasize the importance of a non-zero value of 
Afc/. As discussed by Clarke, Strong and Anderson [ssj. 
the properties of Luttinger liquid, which have been estab- 
lished for a single one-dimensional chain, can extend in 
the case of a quasi-one-dimensional system (spin/charge 
separation, power-law behaviour of correlation func- 
tions); that is, even though the band structure extends 
in the _L dimension (as Akf ^ 0), the system will not 
converge to the two-dimensional Fermi liquid. These au- 
thors claim that _L superconductivity originates from this 
mechanism, which also relates to confinement in the || 
dimension. From this point of view, unconventional su- 
perconductivity and Luttinger liquid concept (in partic- 
ular spin/charge separation) are interplaying; this gives 
an explanation for the possibility of coexistence of SDW 
instabilities and superconductivity. 

In this SC phase, we also observe a quantitative dif- 
ference between the results obtained using a OPI or a 
Wick-ordered scheme. In the first case, the value of Akf 
lies in the interval [25,30], while in the second, it lies in 
[4,5] (see Fig.[7](c) and (d)). 

In the SDW phase, Afc/ ^ as A ^ (see Fig. [7] 
(a) and (b)). This proves that this phase relates to the 
Luttinger solution. Contrary to the SC phase, the band 
structure remains purely one-dimensional. This system, 
however, is different from Luttinger's original one, be- 
cause, in real space, there are still two chains, with non- 
zero hopping in-between. Let us examine in detail the 
behaviour of the scattering susceptibilities. Couplings 
goi and go2 are not diverging (see Fig. [5]), contrary to 
what is observed in the SC phase. The curves of all cou- 
plings are very close to those obtained when Afc/ is kept 
constant. These results are very different from that of 
Fabrizio, who finds the behaviour of a single chain (see 
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Figure 8: Flow of couplings goi and 302 during the RG flow, 
for U = 0.5 and 2t±/Ao = 1.485. 

Ref. [2H and compare with Fig. [5] here). 

One must observe that, as Akf — > 0, the effective range 
of /^ii shrinks, so that, when one approaches the critical 
scale Ac, the K^^ dependence is extremely badly taken 
into account. On could even expect to recover usual RG 
calculations, for which no SDW phase is found; however, 
the observation of this phase has proved surprisingly ro- 
bust; the explanation is probably that, just before Ac, 
the divergence of Xs^ dominates already in such a way 
that it prevents any divergence of Xs*^^''''- Nevertheless, 
this discussion sheds also light about serious numerical 
convergence problems that arise in this region, and have 
required technical answers. 

Let us compare these calculations with experimental 
data. No direct determination of U or t± are avail- 
able, one can only get indirect determinations by match- 
ing experimental and theoretical curves, as it is done 
in Ref. [s^, for Sri4-xCaxCu2404i compounds (with 
X = 12). These authors compare several experimen- 
tal and theoretical spin susceptibility curves (including 
uniform spin susceptibilities) and obtain a best fit for 



U/t^^ ~ 4 and \t±\ ^ t|| (those values correspond here 
to U ^ 1 and 2t± ~ Aq), which corresponds here to 
predictions done with the Wick-ordered scheme; other 
determinations are available, which correspond to the 
OPI scheme. Actually, the determinations are accurate 
within an order of magnitude, thus it is not possible, 
from experimental data, to decide which scheme gives 
the right predictions. However, many theoretical pre- 
dictions, and a few experimental fits, are located in this 
region of parameters, close to point (J7 ~ 1, 2t±/Ao ^ 1) 
in the phase diagram (for instance by a factor 2). In 
particular the boundary, between SC and SDW phases, 
is located in the same region of parameters. Therefore, 
even if we can't discriminate between Wick-ordered and 
OPI schemes, the phase diagrams we have calculated are 
qualitatively in good agreement with experimental obser- 
vations. 

In conclusion, we would like to emphasize that these 
calculations confirm the determination of a pure SDW 
phase using a very simple ladder model, which was far 
from being obvious until now. The SC phase also indi- 
cates a possible coexistence of magnetism and supercon- 
ductivity, as it is indeed observed, both theoretically and 
experimentally. 

We have also established the importance of the choice 
of the RG scheme. Even if this alternative only raises 
quantitative differences, they are not negligible, so this 
has to be carefully taken into account. We hope that in 
further and more precise models, a clear discrimination 
between the two schemes will be possible, and that it 
will confirm our conjecture that the OPI scheme is more 
accurate. 

The authors would like to thank deeply Christoph 
Nickel for his collaboration in the theoretical preparation 
of this paper; without his help, this work would probably 
never have been completed successfully. 



[1] M. Azuma, Z. Hiroi, M. Takano, K. Ishida & Y. Kitaoka, 

Phys. Rev. Lett. 73, 3463 (1994). 
[2] Le Canon, Mater. Res. Bull. 23, 1355 (1988); Sigrist, 

ibidem, 1429. 
[3] Z. Hiroi & M. Takano, Nature 377, 41 (1995). 
[4] H. Sato & M. Naito, Physica C 274, 221 (1997). 
[5] E. Dagotto, Rep. Prog. Phys. 62, 1525 (1999). 
[6] Y. Piskunov et al, EPJ B 13, 417 (2000). 
[7] See Fig. 3 in: J. Hopkinson & K. Le Hur, Phys. Rev. B 

69, 245105 (2004). 
[8] K. Kuroki, T. Kimura & H. Aoki, Phys. Rev. B 54, 

R15641 (1996). 
[9] J. A. Riera, Phys. Rev. B 64, 104520 (2001) 
[10] M. Randeria, N. Trivedi, A. Moreo, R. T. Scalettar, Phys. 

Rev. Lett. 69, 2001 (1992). 
[11] N. Bulut, D. J. Scalapino & S. R. White, Phys. Rev. B 

50, 7215 (1994). 
[12] P. Germain & M. Lagus, EPJ B 8, 497 (1999). 



[13] R. M. Noack, S. R. White & D. J. Scalapino, Phys. Rev. 

Lett. 73, 882 (1994). 
[14] G. Sierra & M. A. Martin-Delgado, Phys. Rev. B 56, 

8774 (1997). 

[15] A. Moreo, S. Haas, A. W. Sandvik & E. Dagotto, Phys. 
Rev. B 51, 12045 (1995). 

[16] K. G. Wilson, Phys. Rev. B 4, 3174 (1971). 

[17] K. Kuroki & H. Aoki, Phys. Rev. Lett. 72, 2947 (1994). 

[18] H. J. Schulz, Phys. Rev. B 53, R2959 (1996). 

[19] T. Giamarchi & H. J. Schulz, J. Phys. (Fr) 49, 819 (1988). 

[20] It is impossible to review the hundreds of articles us- 
ing this model; three pioneer works are P. W. Anderson, 
Science 235, 1196 (1987), J. R. Schrieffer, X.-G. Wen 
& S.-C. Zhang, Phys. Rev. Lett. 60, 944 (1988) and C. 
L. Kane, P. A. Lee & N. Read, Phys. Rev. B 39, 6880 
(1989). 

[21] M. Fabrizio, Phys. Rev. B 48, 15838 (1993). 

[22] G. Abramovici, J. C. Nickel & M. Heritier, Phys. Rev. B 



7 



72, 045120 (2005). 
[23] W. Metzner, C. Castellani & C. Di Castro, Adv. Phys. 
47, 317 (1998). 

[24] C. Honerkamp, M. Salmhofer, N. Furukawa & T. M. Rice, 

Phys. Rev. B 63, 35109 (2001). 
[25] M. Salmhofer, Phys. Lett. B 408, 245 (1997), M. 

Salmhofer, Commun. Math. Phys. 194, 249 (1998). 
[26] J. Solyom, Adv. Phys. 28, 201 (1979). 
[27] N. Dupuis, Eur. Phys. B 3, 315 (1998). 
[28] S. Dusuel & B. Dougot, Phys. Rev. B 67, 205111 (2003). 
[29] D. Rohe & W. Metzner, Phys. Rev. B 71, 115116 (2005). 
[30] A. A. Katanin & A. P. Kampf, Phys. Rev. Let. 93, 106406 

(2004). 

[31] For a pedagogical explanation of the Wick-ordered 
scheme, see C. J. Halboth, Als Ms gedr.-Aachen: Shaker 
1999, ed. Techn. Hochsch. (Aachen), Diss., 1999 (ISBN 
3-8265-6822-2). 

[32] For a pedagogical explanation of the OPI scheme, see J. 
C. Nickel, These de troisieme cycle, de I'Universite Paris 
11 (2004). 

[33] Other choices are possible, see for instance the scheme 



used in: C. Bourbonnais and R. Duprat, Bulletin of the 
American Physical Society, 49 No. 1, 179 (2004). 

[34] In one-loop diagrams, there is only one integration, how- 
ever, the domain of external momenta are different in the 
OPI scheme and in the Wick-ordered one. 

[35] All calculations in Ref. [2^ are done using the OPI 
scheme. They have been also performed using the Wick- 
ordered scheme, but there are no noticeable differences. 

[36] D. Zanchi, Europhys. Lett. 55, 376 (2001), D. Zanchi & 
H. J. Schulz, Phys. Rev. B 61, 13609 (2000). 

[37] It is a fact that all curves are exactly converging to the 
point (2tx/Ao = 0.02, f/ = 0), however, one has to take 
these extrapolations very carefully, because the theoret- 
ical precision of the calculations behaves like Ac/Ao and 
one observes that log(Ac/Ao) > 100 for 2t±/Ao < 0.01. 

[38] D. G. Clarke, S. P. Strong & P. W. Anderson, Phys. Rev. 
Lett. 74, 4499 (1995). 

[39] See last §(before acknowlegments) in: M. Tsuchiizu & Y. 
Suzumura, Phys. Rev. B 72, 075121 (2005). 



